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ON 



We analyse high-field current fluctuations in metallic systems by direct mapping of the Fermi- 
liquid correlations to the semiclassical nonequilibrium state. We give three applications. First, 
for bulk conductors, we show that there is a unique nonequilibrium analogue to the fluctuation- 
dissipation theorem for thermal noise. With it, we calculate suppression of the excess hot-electron 
term by Pauli exclusion. Second, in the degenerate mesoscopic regime, we argue that shot noise and 
thermal noise are incommensurate. They cannot be connected by a smooth, universal interpolation 
formula. This follows from their contrasting responses to Coulomb screening. We propose an 
experiment to test this mismatch. Third, we carry out an exact model calculation of high-field shot 
*^ noise in narrow mesoscopic wires. We show that a distinctive mode of suppression arises from the 

C/3 ' structure of the Boltzmann equation in two and three dimensions. In one dimension such a mode 

D does not exist. 

Ej: 

I. INTRODUCTION 

s . 

' High-field noispj in degenerate conductors still lacks a systematic theoretical description, despite its importance for 



microelectronics. lil Here we advance a practicable theory of nonequilibrium fluctuations in metallic systems, accounting 
for the dominant Fermi-liquid behaviour of the electrons .□ There is no electronic property of a metal near equilibrium 



a 

q ' that is not governed by Pauli exclusion, from the microscopic level to the bulkE In this paper we analyse how 
degeneracy determines nonequilibrium current noise over a wide range of length scales, even at high fields. 

■ Technological developments of late have led to a variety of delicate measurements of transport aad_noise in many 
| different mesostructures.Qllj Alongside the experiments there has been much theoretical activity.ElTEZl A particular 

0/^ ■ topic, still being elaborated, is the behaviour of current fluctuations in mesoscopic conductors. These are typically 
' shorter than the inelastic mean free path but longer than that for elastic processes. They are in the regime of diffusive 

C*") . transport,, where two diverse understandings predominate. One technique phAe quantum-transmissioit-ZLandauer) 

ON ■ methodMl as applied to fluctuations by Lesovik,E3 Beenakker and Biittiker Martin and Landauerjia and many 
, othersO The secondpao«roach is the semiclassical-transport (Boltzmann) methodES associated with NagaevE3 and de 
Jong and Beenakker .EiltZl Both view a mesoscopic wire a& an assembly of individual elastic scatterers in a bath of free 

-» — \ carriers. The formalisms, however, are quite dissimilar .E9 

■ In the Landauer model, multiple scattering preserves the coherence of single-particle propagation; the only way in 
| which current fluctuations can lose correlation strength is by interplay-pf the transmission amplitudes and nonlocal 

I . Pauli blocking between state occupancies in the source and drain leads.Eil In the Boltzmann model, incoherence enters 
'"O ' from the start through the Stosszahlansatz. To study the fluctuations, the stochastic collision term is supplemented 
Ch , with a set of Langevin flux sources whose phenomenology is that of classical shot noise;Ej their self-correlations reflect 
the sporadic timing of elementary encounters between discrete wave packets and scatterers. The uiduced fluctuations 
. . | lose correlation strength when diffusive elastic scattering is locally modified by Pauli exclusion. OE§Ea 

Regardless of their differences, both methods can explain observations such as the threefold suppression of shot 
noise when elastic scattering prevails.cl Because the phase-coherent model is fully quantum mechanical, its fluctuation 
structure is natural, not imposed, and its predictions enjoy a definitive status. On the other hand, semiclassical 
d ' phenomenology is the more natural tool when inelastic (hence irreducibly phase-breaking) collisions are important, 
as in high-field transport. 

Our paper investigates fluctuations beyond the elastic weak-field limit. A mesoscopic sample is easily driven into the 
high-field regime; some tens of millivolts across a length of 100 nm will do it.Ej Nevertheless, although nonequilibrium 
noise is a unique source of dynamical information out of reach to linear-response theory£3 none of the existing models 
has been pushed substantially beyond its low-field perturbative regime. The Landauer and Boltzmann-Langevin 
formalisms each suffers from its own obstacles to addressing strongly nonequilibrium effects, as does the quantum- 
kinetic theory of Altshuler, Levitov, and YakovetsJHS which tries to unify them. There is, therefore, a real need for 
another approach: one that is nonperturbative in the driving field. This need has been reasserted very recently by 
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numerical evidence of a rich structure for shot noise at high fields, even in nondegenerate systems.E3 

We propose a direct andrViejsatile method for treating noise semiclassically from equilibrium up to high fields, build- 
ing on Fermi-liquid theorjoES and a family o£ Green functions for the linearised transport equation. Such functions 
have been studied by Kogan and Sh|Ul' jman,Ej Gantscvich, Gurevich, and Katilius£3 and by Stanton and Wilkins, 
in great detail, for semiconductors] 2 ^P ! Unlike Boltzmann-Langevin, this approach is not limited to the specifically 
Boltzmannian form of the collision integral.E3 For every collisional approximation that can be used in the one-particle 
transport equation, there is a systematic construction for the two-particle Boltzmann-Green functions (BGFs). These 
then generate both the steady-state and transient nonequilibrium fluctuations, respecting the conservation laws as 
well as the one-body collisional structure and bringing great flexibility to noise calculations. An application to the 
noise perfomance of heterojunction transistors is given by Green and ChiversO __ _ 

In Section II we present a general framework for noise in small metallic systems, including Coulomb effectsEj'EScil 
where we identify two distinct mechanisms for screening. Sec. Ill contains the major applications of our theoijy- 
First, we extend the fluctuation-dissipation relation for thermally driven noise to nonequilibrium-|Conductors,E3~EII 
highlighting the role of Pauli exclusion in suppressing hot-electron effects in the bulk noise spectrum.^! We then discuss 
the many-body origin of shot noise and argue that, in their sharply contrasting responses to Coulomb screening, shot 
noise and thermal noise display quite different physics. We propose a simple experimental test of this difference. In 
Sec. IV we make an exact computation of shot noise within the Drude picture of a conducting wire. We include the 
effects of finite wire thickness on carrier motion, a significant source of shot-noise suppression at high currents that is 
unrelated to diffusive elastic scattering and to Coulomb screening. We sum up in Sec. V. 



II. THEORY 



The theoretical discussion is in four parts. We begin by formulating the transport problem as a direct mapping of 
the electron Fermi liquid to its nonequilibrium steady state. Second, we describe the steady-state fluctuations. Third, 
we discuss the dynamic fluctuations and their formal connection with the steady state. Last, we incorporate Coulomb 
screening into the nonequilibrium structure. Our end product is a complete expression for the current autocorrelation, 
which determines the noise. 



A. Transport Model 



The semiclassical Boltzmann transport equation (BTE) for the electron distribution function f a (t) = / s (r,k, t) is 
d 



d e E(r,t) d_ 
U 9k 



fa(t) = - £ [W*> a (l ~ foe') fa ~ W aa ,(l - f a )f a , 



(1) 



Label a = {r, k, s} denotes a point in single-particle phase space, while sub-label s indexes both the discrete subbands 
(or valleys) of a multi-level system and the spin state. The system is acted upon by the total internal field E(r, t). We 
study single-particle scattering, with a rate W aa < = 5(r— r')W ss < (r, k, k') that is local in real space, independent of the 
driving field, and that satisfies detailed balance: W a > a {\ — f^)fa q — W aa '(l ~ fa q )fa^ wnere fa q 1S the equilibrium 
distribution. In a system with v dimensions, we make the following correspondence for the identity operator: 



5rr> 



{fi(r)5 kk ,} «— > S ss ,S(r - r')(27r)"S(k - k'). 



The volume fi(r) of a local cell in real space becomes the measure for spatial integration, while its inverse defines the 
scaling in wave- vector space for the local bands {k, s}. 

The first step is to construct the steady-state solution f a = f a (t — > oo) explicitly from / cq , which satisfies the 
equilibrium, collisionless form of Eq. (m): 



df? eE (r) df? 
Vk s -— ^— = U. 



Or 



dk 



(2) 



The internal field Eo(r) is defined in the absence of a driving field. The quantities / oq and Eo are linked self-consistently 
by the usual constitutive relations, the first being the Poisson equation 



JUe,, = -47re| „ (r 



(3) 
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in terms of the dielectric constant e(r), the electron density (/ oq (r)) = fKr) 1 J^ k s /° q , and the positive background 
density n + (r) , which will remain unchanged throughout our calculations.^ Normalisation to the total particle number 
is J2 r ^( r )(/° q ( r )) = N. The second relation is the form of the equilibrium function itself, 



feq 
J a 



1 + exp 



k B T 



(4) 



in which the local conduction-band energy e a = e s (k; r) can have band parameters that depend on position. The 
locally defined Fermi level <f> a — /i — Vo(r) is the difference of the global chemical potential /i and the electrostatic 
potential Vb(r), whose gradient is eEo(r). 

Define the difference function g a — f a — f° q . From each side of Eq. (|l|) in the steady state, subtract its equilibrium 
counterpart. We obtain 



dg a eE(r) 8g a e(E - E ) df^ ^ 



ft 9k 

a' 

- W aa ,)(fjg a + g a >f? + g a g a > 



(5) 



The solutions to Eqs. (0) and (g) are determined by the asymptotic conditions in the source and drain reservoirs, be 
it at equilibrium or with an external electromotive force. Our active region includes the carriers within the source and 
drain terminals out to several screening lengths, so that local fields are negligible at the interfaces with the reservoirs. 
In practice we assume that all fields are shorted out so that E(r) = Eg(r) = beyond the boundaries. Gauss's 
theorem implies that this bounded system remains neutral overall: 

5>(r)( 5 (r))=£ 5Q = 0. (6) 

r a. 

We put Eq. (|^) into integro-differential form, with an inhomogeneous term explicitly dependent on f cq : 

E S ^ A /]-'5a' = ^-^ + E W -'^ (7) 

ol' a' 

where E(r) = E(r) — E (r) is the local field induced by the external electromotive potential and B[W A f] is the 
linearised Boltzmann operator 



B[W A f] aa , = S aa 



d eE(r') 3 ^. Jir nrA 



dr' h 

W aa ,+W A a ,f a , 



(8) 



and W A a , = W aa i — W a > a . Note that W A = if the scattering is elastic. To represent the physical solution, g must 
vanish with E in the equilibrium limit. This is guaranteed by the Poisson equation 



d_ 
dr 



•eE 



-Aire 



((/(r)) - (rV))) =-47re( 5 (r)>. 



(9) 



B. Boltzmann-Green Functions 



To calculate the adiabatic response of the system about its nonequilibrium steady state we introduce the Boltzmann- 
Green functionEj 



def Sga 



u a« — x rcq j \ ±KJ ) 



with a global constraint following directly from Eq. 
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y ' G aa > — 0. 

a. 

The equation of motion for G is derived by taking variations on both sides of Eq. (||): 



(11) 



Y / B{W A fUG l3a , 



eE(r') d ^ TirA 

— h > Wn , 

h dk> /3q 



W A a ,g 



(12) 



The variation is restricted by excluding the reaction of the local fields Eo(r) and E(r). This means that G is a response 
function free of Coulomb screening. Here we treat the electrons as an effectively neutral Fermi liquid; in Section IID 
we give the complete fluctuation structure, with Coulomb effects. 

All of the steady-state fluctuation properties induced by the thermal background are specified in terms of G and 
the equilibrium fluctuation A/ cq . The equilibrium fluctuation is the proper particle-particle correlation in the static 
long-wavelength limit, normalised by the thermal energy ksT . In the Lindhard approximation!] this is 



A/^ = fc fl T^ = /Hl-/* q ). 
o<p a 



(13) 



When there are strong exchange-correlation interactions, this two-body expression is renormalised by a factor depen- 
dent on the Landau quasiparticle parameters.00 In this work we consider free electrons only. 

Define the two-particle fluctuation function A/^ a , = (S aa i + G QQ <)A/°?. The steady-state distribution of the 
particle-number fluctuation is the sum of all of the two-body terms: 



Af a = J2 A/5 = Af? + J2 G aa <AfJ. 



(14) 



The noncquilibrium fluctuation A/, manifestly a linear functional of its equilibrium Fcrmi-Dirac form, is the exact 
solution to the linearised Boltzmann equation: 



Y,B[W A fUAf 



0. 



(15) 



Charge neutrality implies that the total fluctuation strength over the sample, AiV = ^Z r f2(r)(A/(r)), is conserved. 
This constrains both steady-state and dynamical fluctuations. ■— ■ 
Calculation of the dynamic response requires the time-dependent Boltzmann-Green functionlij 



Ofa'{t) 



(16) 



with initial value R aa > (0) = S aa i . As with G, the variation is restricted. The linearised BTE satisfied by R(t — t') is 
derived from Eq. ([!]) and takes the form 



[ 5 "0§- t + B[W A f] a ^Rp a ,{t - f ) = S(t - t')5 aa ,. 
Summation over a on both sides of this equation leads to conservation of normalisation:0 

Y,R aa -{t-t') = e(t-t>). 



(17) 



(18) 



The time-dependent BGF is a two-point correlation. It tracks the history of a free electron in state a' added to the 
system at time t'\ the probability of finding the electron in state a at time t is just R aa '(t — t'). In the long-time 



limit Eq. 
equation. 



19 



17p goes to its steady-state form, and so R aa '{t oo) oc A/ Q , the solution to the steady-state linearised 
Together with Eq. (Ill) this gives the identity 



R aa , (t -> oo) 



Afa 

AN' 



(19) 
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All of the time-dependent fluctuation properties induced by the thermal background are specified in—terms of 
R and the steady-state nonequilibrium fluctuation A/. From the dynamical two-particle distribution,E£l that is 

(2) 

A/^(f) = R aol '(t)Af a ', one can construct the lowest-order moment 

A/ tt (i) = (20) 

a' 

in analogy with Eq. ()l4|). Equation ([l?]) has an adjoint,@ with whose help one can show that Af a (t) = Af a for t > 0. 
Thus the intrinsic time dependence of Af^it) is not revealed through this quantity. [Note: a remark in Green, Ref. 
|39| , that Af(t) is inherently time-dependent, is true only for collision-time approximations.] Eq. ( fL8[ ) implies that the 
total fluctuation strength is constant: ^ r Q(r)(A/(r, i)) = AN. 



C. Dynamic Correlations 

We move to the frequency domain. An important outcome of this analysis is the extension of the fluctuation- 
dissipation (FD) relation to the nonequilibrium regime. This requires expressing both the difference function g and 
the adiabatic Boltzmann-Green function G in terms of the dynamical response. The Fourier transform R(w) of the 
time-dependent BGF satisfies 

E {B[W A f] af} - iu5 a0 }R 0a ,(u) = <W, (21) 


making R(w) the resolvent for the linearised operator of Eq. (||). The global condition on R(w) from Eq. ( p"8| ) is 

^R^H = - 1 rj^0+. (22) 

a \ i/ 

At first sight this fails to match the corresponding condition on the adiabatic Boltzmann-Green function, Eq. (|Tl|). 
To determine the solution of Eq. (|l2|) for G in terms of the resolvent, we follow Kogan and Shul'martJ and introduce 
the intrinsically correlated part of R(cj), namely 

C aa , H = R QQ , H + 1 ^ . (23) 
i(uj + it]) AN 

This correlated propagator satisfies a pair of identities in the frequency domainjll the transform of the relation 
Af{t) = 9(t)Af leads to 

^C m -HA/ a ,=0 (24a) 

a' 

while Eq. ( |22| ) leads to 

E<wH=0. (24b) 

a 

The second of these corresponds to the constraint on G. Like R(w), the correlated BGF is analytic in the upper 
half-plane lm{o;} > 0, and satisfies the Kramers-Kronig dispersion relations. Unlike R(w), however, C(oj) is regular 
for uj — > 0. 

We now obtain g and G in terms of the correlated dynamical propagator. Consider the equation 

E { B {W A fW - io;daa'}ga'M = + E^ W -'^'> (25) 

a' a' 

with solution 

&» - E C -'^)^- i '^T + E Caa'Mffa'W.W (26) 
a' a'/3 
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The uncorrelated component of R(w) does not contribute to the right-hand side of this equation; in the first term it 
results in a decoupling of the summation over a! ', yielding zero because df^/dk.' is odd in k'. In the second term, 
decoupling means that the double summation over a' and f3 vanishes by antisymmetry. In the static limit Eq. ( |25| ) 
becomes the inhomogeneous equation (Q), and moreover g(0) satisfies Eq. (||), the sum rule for g. Therefore 

pV(r') 8f cq 

9a =E C -'(°)^- i -^T+E C -'(°)^'^V/3- (27) 

a' a'f3 

This identity is central to the FD relation. 

In models with symmetric scattering W is zero and the adiabatic BGF assumes a simple form on varying both 
sides of Eq. @: 

G aa ' = C QCe '(0) — - — "g^i- (28) 

More generally, an analysis similar to that for g(w) can be used directly for the adiabatic propagator. Introduce the 
operator G(ui), defined to satisfy the dynamic extension of Eq. (Oh, 



{B[W A fU - ico8 afi }Gp a ,(u;) = S a 

13 



eE(r') d yr^ TIrA 



W£ a ,g a . (29) 



G aa >(u) = C aa ,(u;)^p-~ - (Caa>(u) - C a p{u)j)w£ g p . (30) 



This has the solution 

dk> 

p 

In the first term on the right-hand side, the uncorrelated component of R(w) makes no contribution after decoupling 
because the physical distributions F a on which G(w) operates vanish sufficiently fast that ^2 k dF a /d\t = 0. In the 
second right-hand term the uncorrelated parts of K aa ' and K a/ 3 cancel directly. We conclude as before that 

G aa , = C -'(°)^y^-^7 - E ( C ««'(°) - C<*p(0j)w£f,gf,. (31) 



This is a crucial result because it shows (a) that the adiabatic structure of the steady state, through G, is of one 
piece with the dynamics, and (b) that the nonequilibrium correlations originate manifestly from the equilibrium state, 
through GA/ cq . We have thus proved that the kinetic BGF analysis is formally self-sufficient once its boundary 
conditions are given. This fact is embodied in the frequency sum rule 

™ = - f-HG^M}. (32) 

The vehicle for the physics of current noise is the velocity autocorrelation. It is a two-point distribution in real 
space, built on the correlated part of the two-particle fluctuation Af!^, = R aa > Af a > and taking the formEl 

«vv'Af( 2 )(r,r»)): d ^ f -L _J_ £ £ v ks Re{C QQ < (u;)}v w Af a , (33) 

^ ' ^ ' k,s k',s' 

The nonlocal velocity autocorrelation provides the direct basis for shot-noise calculations when the distance |r — r'| 
becomes comparable to the mean free path. The local function derived from it, 

Sfaw) = e 2 ^fi(r')(((E(r).v)(E(r')-v')Af( 2 )(r,r';^)))' c , (34) 

has a macroscopic reach since in effect it samples fluctuations over the bulk. It is closely related to the current-noise 
spectral density and satisfies the nonequilibrium FD relation discussed in Sec. III. 



6 



D. Coulomb Effects 



We generate the Boltzmann- Green functions in the presence of induced fluctuations of the electric fields. Variations 
are now unrestricted. The resulting description of screening effects extends the physics of classical space-charge 
suppression.^ We consider samples with a fixed dielectric constant, and likewise for the reservoirs. 

The equation for the screened equilibrium fluctuation A/ eq = kBTSf eq /5pL is obtained by operating on Eq. (§) 
satisfied by the equilibrium distribution. We have 

dAf^ eE (r) gAjg ( e ^ 5E (r) ~ , eq \ df? 

^ s ~ r n diT= [h^f^r A - L ')-^k- (35) 

Detailed balance keeps the equation collisionless while the Poisson equation (||) implies that, within the system 
boundaries, the variation of eE with respect to / eq is the Coulomb force for an electron, 

e^fp? = ~eEc(r - r') = ^V c (r - r') (36) 

where Vc(t) = e 2 /e|r| is the Coulomb potential. As a result Eq. ( |35| ) becomes 

a' 

Viewed as a variant of the equilibrium BTE, Eq. d35|') is inhomogeneous. Its solution includes a term proportional to 
the homogeneous solution, which in this case is the bare fluctuation A/ cq . Let 7c be the proportionality constant. 
Then 

A/° q = >y c Af? -JE C^,(0)^-E c (r' - r")A/^„ (37) 

a 1 ol" 

in which C oq is the correlated part of the resolvent for the equilibrium state. The integral on the right-hand side of 
Eq. ( p7| ) has a structure similar to Eq. (27), in that the uncorrelated part of the resolvent gives no contribution after 
decoupling of the intermediate wave- vector sums. 



1. Thomas-Fermi Screening 

The constant 7c is sensitive to the physics of charge transfer between sample and reservoirs. Recall that the 
fluctuation A/° q = ksTdf^ 1 / d(j) a is a measure of the electrons' response, as a Fermi liquid, to a change in the 
effective Fermi level <p a = /1 — Vo(r). The latter is the net contribution from kinematics alone to the cost of adding 
an electron locally to the system; the electrostatic energy Vb(r) is excluded from the Fermi- liquid accounting. 

When the Coulomb fields are frozen, as in the restricted analysis, the Fermi-level variation is that of the global 
chemical potential: 5<p a — 8fi. In the full Coulomb problem, we must offset the energy cost of charge transfer from 
reservoir to sample. The Coulomb energy needed to add an electron to the conductor is 

a 

A corresponding term u r characterises the reservoirs. However, u r cannot be probed directly; its effects are absorbed 
within the operational definition of the chemical potential. This means that u r = identically, and that u c thus 
represents the net work to move an electron from reservoir to sample. It is the conduction-electron contribution 
to the contact potential.^ (By contrast, the core-electron contribution determines the offset in the band bottom 
£ s (k=0; r). This is independent of y, and does not appear explicitly in the variational analysis.) 

It follows that the portion of the chemical potential sustaining the Fermi liquid in the conductor is \i — u c . Free 
variation of the global parameter /i generates the coefficient 
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Application to Eq. ( |37| ) of Eq. (24t), that is the sum rule ^2 a C aa > — 0, establishes the normalisation 

^A/^ = 7c AiV. (40) 

a 

If sample and reservoir have matching electronic properties, then u c = u r = and 7c = 1. This is the norm for noise 
measurements in metallic wires. If, on the other hand, the sample differs substantially from the reservoir in metallic 
structure, then from Eq. ( |38|) and the form for Vq(t) as the solution to the Poisson equation (||) one obtains 

™ = 1 - jf E f E Vc{r - r<)C + V (r) - |) ^J, (41) 

a \ q' / 

an instance of suppression by self-consistent Thomas-Fermi screening. 

There is strong indirect evidence that this mechanism is the major determinant of low-noise performance in het- 
erojunction field-effect devices£j At mesoscopic scales, one can anticipate a wide variety of interfacial screening 
behaviours for the fluctuations. Potentially interesting is the case of metal-semiconductor-metal structures. See 
Appendix Eq. JaTo| ), 

Note that 7c enters only at the two-body level. It cannot renormalise the one-body distribution g, or any averages 
constructed with g. This includes transport coefficients such as the mobility. 

Equation ( |37| ) is solved by introducing a Coulomb screening operator r eq (0), whose inverse is 

eq 



r-fl))- 1 ) = <W + C ^(0) ag--Ec(r, - r'). (42a) 

This yields 

A/« q = EC(0)(7cA/:?), (42b) 

a' 

analogous to the Lindhard screening theory of the electron gas in the static limit .i Taken together with Eq. (|4l|), it 
allows a closed-form solution for 7c. 



2. Collision-Mediated Screening 
Away from equilibrium we define the Coulomb screening operator through its inverse 

(r(c)- 1 )^ ^ 6 aa , + C a/3 ( W )^-E c (r /3 - r'). (43) 

While 7c is collisionless, the operator T captures the dynamics of interaction between scattering and screening, an 
exclusively nonequilibrium process. A few elementary results for both collisional and Thomas-Fermi screening are 
discussed in Appendix An important property of the collision-mediated screening operator, following from Eq. 
«,is 



E r -'H = L ( 44 ) 

a 

The significance of T first becomes evident in obtaining the screened adiabatic propagator G aa i = Sg a /Sf^}, whose 
unrestricted BTE [cf Eq. (fil) for G] is 



E B{W-f U G Pa/ = E B[W*fUG 0a/ + J (E ) ^ 

e<5E (r) df? 



h 5f3 dk 

Poisson's equation (Q) once again determines the variation of E with respect to / as 



(45) 
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3E(r) 
~5f7 



= -E c (r-r'). 



(46) 



The solution of Eq. (45) is a two-pass process, in which one first resolves the Boltzmann operator B[W A f] and then 
invokes T to rationalise the transformed equation: 



G aa ' = G aa > — — E C Q| a(0)^r-^--Ec(r/3 — 1*3') (<V a' + Gp' a 



^C a ,(0)^.E c (r,-r') 



E r ^(0) Up* ~ |E C ^'(°)^- E c(^ -r')j 



We have also used Sfp'/Sf™? — Sp/ a i + Gpi a i. The screened steady-state fluctuation is now 



(47) 



(48) 



The main outcome of the structure of G is the invariance of the fluctuation strength over the sample; this follows 
from ^2 a G aa ' = 0. Thus 



E A /« = "fcAN ee AN. 



(49) 



3. Dynamics 

We now examine the screened dynamics. In the time domain the screened resolvent R{t) has the same formal 
definition, Eq. (|l6|), as its restricted analogue. The unrestricted equation of motion, Fourier transformed, is 

V {B[W A f} a/3 - kjS aP }R Pa f(u) = 5 aa , - J V ?j^-Ec(r - r^R^H, (50) 

with solution 

R QQ / (w) = 2J T a p(aj)Rf3 a > (u>). (51) 



This resolvent obeys identities analogous to those for R(w), namely Eq. (|22|) and, in the time domain, R aa > (t — + oo) = 



Af a / AN. Together with Eq. fl49|), equality of the residues at lo — on each side of Eq. (51) implies the relation 

Afa = 7c E ^ aa ' (°) ' (52) 



equivalent to Eq. (48) by the properties of G and .BpW^ 4 /]. 
The correlated propagator, with screening, is 



C aa ' (<jj) = R QQ ' (<jj) + —, — ^——^h-. (53) 
liuj + ir]) AN 



Its structure follows from combining Eqs. (p3J), (|5l|), and (b2p for 



|9 /? 
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Applying the screened form of Eq. (24a), that is J2B/C/3p>(w)Afpi = 0, to Eq. ([34|) generates 



00' 

Fed back into Eq. (|54|), this produces 



C Q q'(cj) = ^r Q/3 (cj)C0/5'(cj) ( 5/3> a > - ) • ( 55 ) 



00' 



4- Current- Current Correlation 

All the components are in place to construct the velocity autocorrelation function in the presence of screening. At 
zero frequency this is 

<(vv'Af( 2 )(r,r';0)))' c d = f J-^^]T ]T v ks C aa , (0)v kv A/ a „ (56) 

At finite frequency we must add the displacement-current contribution to the fluctuations. The velocity is replaced 
with the nonlocal operator 

u ks (ry^)^v ks -^E c (r-r''), (57) 

which requires two intermediate spatial sums to be incorporated within the expectation ((Re{uAf( 2 )u'*}))' c . For u> = 
this recovers Eq. (|56|), a more complex expression than its bare counterpart Pfli--|(^3|)- ln practice, collisional Coulomb 
effects are dominant in mesoscopic and in strongly inhomogeneous systems EU'LJ 

III. APPLICATIONS 

The first of our applications connects thermal fluctuations and dissipation in the bulk nonequilibrium context. 
Little is known of the effects of degeneracy on noise beyond the linear limit ,□ and we analyse them here. In our second 
application we investigate the many-body nature of mesoscopic shot noise. For degenerate electrons we show that 
thermal and shot noise have very different physical properties not easily subsumed under a single formula.liil 

A. Nonequilibrium Fluctuation-Dissipation Relation 

The fluctuation-dissipation relation near equilibrium connects the spectral density of the thermal fluctuations to 
the dissipative effects of the steady current in the system. However, dissipation by itself does not exhaust the physics 
of this sum rule. iTJiere are nonlinear terms, negligible in linear response, that dominate the high-field behaviour of 
the current noiseOLy We calculate these contributions. Since the relation is macroscopic, to lowest order we omit 
Coulomb screening effects; these are weak in the bulk metallic limit (see for example Appendix 

The resolvent property of R(w) provides a formal connection between the steady-state (one-body) solution g and 
the dynamical (twp-body) fluctuation Af at the semiclassical level. Taken to its equilibrium limit this becomes the 
familiar theorem.0 The connection is made in two steps. Consider the kinematic identity 

df cq h 

-k = -ur^i- (58) 

and apply it to the leading term on the right-hand side of Eq. (^7|) . The result is 

e 



9o 



k B T 



J2C aa >{0)(E-v) a ,Af°? + h a , (59) 
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in which h a = ^2 a ipC aa ' {fyga'W^iggp. Evaluating the current density according to J(r) = —e(vg), the power density 
-P(r) = E(r)«J(r) for Joule heating can be written as 

P(r) = ^ f ^y^^(E.v) a C a ^(Q)(E.v) a ,A/:?-e(E.v/ l ). (60) 

k,s a' 

In the second step we take the one-point spectral function Sf in the static limit, substituting for A/ from Eq. ( |l4| ) 
in the right-hand side of Eq. (|34| ) to give 

S /( r '°) = ^EEE( i ' v )«C oa '(0)(E.v) a ,A/:? + 5 9 (r,0), (61) 



where S g (r, 0) is generated by replacing A/ with Ag = GAf cq in Eq. (p3[), and subsequently in Eq. (p4[). Direct 
comparison of Eqs. ( |60|) and (|6l]) leads to 

^.P (r) + e(g . vft) + ^2). 

This is the nonequilibrium FD relation. 

The standard linear-response result follows. The term in h on the right-hand side varies as Eg 2 , while the final 
term varies as E 2 Ag; therefore both of these contributions are of order E 3 . Suppose that the system is homogeneous 
and that E = E acts along the x-axis: then division by E 2 on both sides of Eq. (p2[) gives 



1 S,(t,0) . \J B \_ at (63) 



E 2 k B T E 



where o is the low-field conductivity. Eq. ( |63| ) is the near-equilibrium statement. 

The purely nonequilibrium structures beyond P(r) can be expanded similarly to it. We discuss the symmetric- 
scattering case, for which there is no contribution e(E-v/i). Within S g we apply the formula for the adiabatic 
Boltzmann-Green function, Eq. (|28|), to express Ag in terms of the correlated propagator C. This produces two 
equivalent closed forms for the higher-order correlation: 

e 2 ~ ~ pPfr'l f)Af cq 

Ss(r > 0) = ^EE( e ' v )«^( )( &v )^E c ^(°)t^'if' (64a) 

^ ' k,s a' 



S g (r,0) 



k B TQ(r 



• ^ ^(E.v) Q (C(0)E.v) QQ ,(l - 20)AC. 



(64b) 



Equation (64b) follows from (64a) after using Eq. ( p8[ ) to express <9A/ cq /<9k in terms of / eq and A/ oq , and absorbing 
an inner sum into (CE-v) 2 . 

The term above is markedly different from the rate of energy loss P(r) from Joule heating. By contrast, S g (r, 0) 
relates directly to nonequilibrium broadening of the fluctuations, due to the kinetic energy gained during ballistic 
motion; the extent of the broadening is dynamically constrained by dissipation. The impact of this term on current 
noise is felt only for significant departures from equilibrium. 

In a degenerate system ther e is a n additional, purely kinematic, constraint on field-driven broadening, seen directly 



in the factor (1 — 2/ cq ) of Eq. (64b). This inhibits the contribution of S g relative to the corresponding classical result, 



in which the factor is unity. Suppression of electron heating by Fermi-Dirac statistics reflects the large energy cost of 
displacing electrons deep inside the Fermi sea. 

To hiffihJiftht the difference between dissipative and hot-electron terms it is instructive to revisit a simple 
example J^J'cll the uniform electron gas in the constant collision time (Drude) approximation, subject to a field 
E = — E~k. Expressions for the power density P and hot-electron component S g are derived in Appendix M. The 
thermally driven current-current spectral density, over a uniform sample of length L x and total volume 17, isS 



S(£,u,)=4^Q(r)£Q(r') 



ev x 
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= 4 



LIE 2 ' 



(65) 



Introducing the conductance Q = QP / L 2 E 2 , the static limit of the spectrum is determined by Eq. (p2[) 



S{E,0)=4Gk B T 



Sg(0) 



PknT 



<±Qk B T 



An ( m*njE 



*,,2 p2\ i 



k B T 



(66) 



We have substituted for P and S g from Eqs. (B7) and (|BS|). The electronic density is n while An = AN/Q is the 
number-fluctuation density. The effective electron mass is m* and \i e is the mobility. 

The term Sg/PksT on the right-hand side of Eq. (66) is a relative measure of the hot-electron contribution to 
the noise. The inhibiting effect of degeneracy, through An/n, is greatest at low temperature; in terms of the Fermi 
energy £p cx n 2 l v we have 



An ksT dn 
n 



n 



dei 



vksT 
2e F 



(67) 



When ef <C ksT the ratio An/n is unity; the hot-electron term is that of a classical electron gas (low density, high 
temperature), whose high-field behaviour is S = 4Qm* /j, 2 E 2 independently of T. On the other hand, when ksT -C £f, 
the system is strongly degenerate and Eq. ( |66| ) with Eq. ( |67| ) yields 



v ( m*a 2 E 2 
2 I E F 



(68) 



The thermal fluctuation spectrum necessarily vanishes with temperature, but its ratio with the Johnson-Nyquist 
spectral density ^QksT continues to exhibit a hot-electron excess, now scaled by the dominant energy £p. Figure [l] 
illustrates the behaviour of the spectral ratio for a two-dimensional electron gas, as a function of the applied field as 
T is taken from the degenerate limit to above the Fermi temperature Tp = £f / ks- We see the gradual trend towards 
the classical form of Eq. (|66| ) with rising temperature. ■ 

Equation ( p8[ ) may be compared with a perturbative estimate by LandauerEil for the degenerate limit, in which 
the analogous hot-electron contribution is (JCZ/fesT) 2 where SU ~ m*n e Evp is a characteristic energy gain and vp 
is the Fermi velocity. Taken at face value, this suggests that hot-electron effects in the low-T regime can be further 
enhanced by cooling. A series expansion in powers of E does not take into account non-analyticity of the Boltzmann 
solutions in the approach to equilibrium j!3 see also Eqs. (C£) and (C14) of our Appendix Non-analyticity of the 
distribution function /k precludes the reliable calculation of moment averages by expanding away from equilibrium. 

The relevance of non-analyticity to transport physics has been questioned by Kubo, Toda, and Hashitsume.li3 They 
ascribe its appearance to the simplistic treatment of real collision processes by the Drude approximation, despite strong 
evidence by Bakshi and Grossc3 that non-analyticity is generic to Boltzmann solutions. Even in the Drude model, the 
nonperturbative solution produces a physically coherent account of the temperature dependence of noncquilibrium 
fluctuations, while finite-order response theory does not. Such clear qualitative differences between perturbative and 
nonperturbative predictions should be detectable in the nonequihbrpim_noise. 

There exist several alternative generalisations of the FD relation. E£rE3 We mention the best known, which defines 
the nonequilibrium noise temperature T n and is pivotal to the interpretation of device-noise data.t3 This effective 
Nyquist temperature is obtained, for a nonlinear operating point, by normalising S with the differential conductance 
Q{E) = (efi/ 'L 2 )d(v x g) / 'dE such that T n (E) — S(E,0) / <iQ(E)k B , corresponding to the output of small-signal noise 
measurements. Our Eqs. (|60| ) - ( |6~I| ) provide a microscopic framework for computing S in a wide class of degenerate 
systems. Since Q{E) is also calculable, this yields T n . 



B. Shot Noise 



Carrier fluctuations manifest as shot noise when they are induced by random changes in the discrete flux at the 
terminals, rather than by thermal agitation distributed through the body of the conductor. Consider an open segment 
of electron gas between macroscopic leads. For this segment we add up the transient, time-of-flight correlations between 
the current at the source boundary, x = x%, and that at the drain boundary, x — x-±. 

The total shot noise measured across the boundaries is the resultant of two components. One component represents 
the response at the drain terminal to the random entry of electrons from the source reservoir, while the other represents 
the response at the source terminal to the random exit of electrons out to the drain reservoir. Thus 
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Ssh(\x2 - xi\) = S^(x2,xi) - S^(x 1 ,x 2 ) (69) 

where the unidirectional term S'z(xj, X4) correlates the induced flux at Xj with the random inducing flux at Xi, and 
has a structure determined as follows. To begin with, note that the correlated two-particle fluctuation at time t, 
following a spontaneous change SN S » in the population of spin subband s", is 

6 fa 



Coulomb screening is fully incorporated. When a particle is added at x± we have 5N S » = +1; when a particle is 
removed at X2, then SN S " = — 1. The sign of 8N S » determines the sign of the corresponding unidirectional term in 
Eq. ©. 

We next observe that, for each of the active carriers in the segment's population N = J2 S " N s ", the time of arrival 
at the source boundary is uncorrelated with all other arrival times. (In the same way, departure times at the drain 
are mutually uncorrelated.) It follows that is an incoherent sum of transients: 



S^x^^nJ e e / 

s" \a[x—Xj] a'[x' —x{\ ' 



dt (v x )^C aa ,(t)(v x )v-^\5N s . 

oN s t> 



2p2 XN J dVr ^ x ~ x i) J rf^'^'-^)^K<Af( 2 )(r,r';0)))' c , (70) 



where we have used Eq. (M9) with the spin trace 



V" at i^L - - X- Sfa ' SN NAfa 



6N S » 2 *-f SN SN S „ AN 

a s 

Our many-body construction of shot noise rests on two assumptions. The first is ergodicity, after the original 
argument of Schottky: a typical carrier in the segment must enter through the source (cathode) and, eventually, leave 
through the drain (anode). There is no temporal correlation among individual transits. The second assumption is 
that each carrier in the ensemble has distinct roles as both agent and spectator: it generates shot noise, and it is also 
part of the many-body response making up the shot noise. These dual roles are statistically independent. .— ■ 

The phenomenological content of Eqs. (|6^) and ( f70| ) is the same as for the Boltzmann-Langevin formalismEj except 
that there is no longer any need for commitment to a specific collisional form (other than expecting it to satisfy 
conservation). Notice too that the total shot noise vanishes identically at equilibrium because detailed balance renders 
the equation of motion for the resolvent R eq (t) collisionless, and hence self-adjoint. Self-adjointness is preserved for 
R cq (t) since Coulomb forces are conservative. One can then show that 

<KRe{C^,MKArV)}>' = «A/ cc »(r)^Re{C a q a ,(-c)}<))'. 



This produces exact cancellation between the right-hand terms of Eq. Q6S^ 

Significantly, the Thomas-Fermi screening coefficient 7c dividing the flux autocorrelation in the second line of Eq. 
([70j), is cancelled exactly by its presence within the autocorrelation via Eqs. (5^) and (|5^). This means that the only 
type of Coulomb screening affecting the shot noise is collision-mediated. Indeed, any homogeneous renormalisation 
of the equilibrium fluctuations leaves the ratio A/ / AN untouched. Consequently such a rescaling has absolutely no 
influence on the shot noise. By comparison, the effect on the free-electron Johnson noise can be dramaticEil since it 
is proportional to 7c- 

The contrast between their Coulomb responses is one demonstration that thermal noise and shot noise are in fact 



distinct many-body phenomena. Although they share a common microscopic structure in Eq. (56), in thermodynamic 
terms one is an extensive continuum quantity driven by fluctuations of the kinetic energy, while the other is short- 
ranged and corpuscular, driven by fluctuations of the local particle number. For strongly degenerate systems the 
two are disproportionate because of the scale difference N/AN, which becomes unity only in the classical limit. It 
follows that, unlike the perturbative treatments, til this formalismj-will not admit a universal interpolation formula 
giving thermal noise in one regime and (true) shot noise elsewherea 

We propose a simple experimental test of incommensurability, applicable at any current. In a point-contact o 
striction defined on a two-dimensional electron gas at a heterojunction, thermal and shot noise are both measurable 
Thermal noise, by its scaling with 7c, depends strongly on electron densityEd Shot noise does not share this depen- 
dence. If the density in the channel is changed, for example by back-gate biasing, the thermal noise should vary 
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strongly with the bias voltage. By contrast, the shot noise should have none of this variation since it is immune to 
self-screening of the carrier fluctuations in the quantum well. 

Our direct concern is the high-field limit, where inelastic collisions rule. For an initial look at high-field shot noise 
we explore the Drude model of a uniform wire, emphasising that its low-field behaviour, though revealing, does not 
address the elastically dominated diffusive regimeJld We take the case where the segment and its leads are physically 
identical.Ej This is an idealised example; experimentally there must be some differentiation between sample and 
reservoirs, expressible in known boundary conditions. 

The region has length I = X2 — x\ <C L x , where L x is the length of the complete assembly, segment plus leads. Since 
at least one of the dimensions may approach the mean free path, the BGFs for this problem have the short-range 
spatial structure detailed in Appendix £J and modified by collisional Coulomb effects. While we do not ana lyse the 
latter in computational detail here, we propose a useful approximation based on the Ansatz of Eq. (A4) for the 
operator T applied to Eqs. ( |52| ) and (jolf). In a conductor of diameter much greater than the mean free path, we take 
the bulk Fourier coefficients C^ 6 ^ (q, 0) and 7 C oii( q, 0) , respectively for the correlated propagator of Eq. ( p!2 ) and the 
collision- mediated Coulomb suppression of Eq. flA2|). The unidirectional shot noise of Eq. (Foh then becomes 



X7coii(q>o) 



/d v q 
^— exp{i[q x (xj - Xi ) + q_L-(r ± - r' ± )]} 



2d v k 



d v k' 



>*)kC^,(q,0)K) k ,A/k' 



(271-)" J (2tt) 



(71) 



where, for any wave vector u, we write its transverse component as ui. We have dropped a contribution ~ Jcou 
coming from the second term on the right-hand side of Eq. ( |55| ) for C. In a full study of nonequilibrium Coulomb 
processes within shot noise, 7 co n is clearly central. Note that there are no transverse terms in one dimension (ID), 
nor is it strictly possible to discuss semiclassical Coulomb effects in ID. 



IV. CALCULATIONS 



In this Section we present calculations for our inelastic model. We omit Coulomb screening, letting 7 co n — > 1 in Eq. 
(fnj). To build up a detailed picture we start with shot noise in a one- dimensional wire. 



A. One Dimension 



For a segment m uch s horte r tha n the total system size, we can simplify the calculation by setting C^ h \q, 0) ~ 
C(°)(g, 0); see Eqs. ( |Cll| ) and flC12j ). The omitted term is proportional to the finite resolution function <pi(q; \L X ) ~ 
5{q) I L x . Since, over the segment, most of the structure involves ql£l, the approximation results in a negligible error 
of order l/L x <C 1. 

Proceeding from Eq. ( |7l|) with v = 1 and using Eq. (|C9b| ) for C<°) we obtain 



lAr 



dq 
I — exp 
2tt 



= 2- 



n 



dk 



iqUl 



f 7rA?i / k, 



k'dk'^6{k - k')e- (k ~ k ' )/kd Af k , 
kd 

1,7 (fc 2 -fc' 2 ) 



2m*kd 



^ J kdke- k / kd J k'dk'e k ' /kd Af k ,s(z - (k 2 - fc' 2 )/^). 



(72) 



We have introduced £ = (xj — Xi)/l for i, j = 1,2 and the wave n umber pd defined by p 2 , = 2m*lkd/TiT — 2m*eV/h 2 . 
Next, use the expression for Af in terms of A/ eq ; see Eq. (C14a). Rearranging the order of integration, we get 



5 a(0 = 2 



ne 2 r f h 2 p 2 d 



m*l V 2m 
= 2el 



dk" Af$ ek „ /kd 



^- e - k ' kd f lkl 2k'dk'5{k' 2 + ip 2 d - k 2 ) 

d J\k"\ 



°°dk"Af$^ k 



(i 



7T 



A? 



"A, £ ^ e -^e{IC)6(k - ^k" 2 +p 2 ) 
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+ e -k"/k d kdk c ~k/k d 
-k" k \ 



(73) 



where / = {ne 2 T /m*l)V is the current through the wire. We now go to the degenerate limit, replacing the equilibrium 
fluctuation according to A/£ q = (m* ksT / Ti 2 kp)8(\k\ — kp) . The number-fluctuation density is An = 2m*kBT jfi'i^kp. 
A series of manipulations leads to the total shot noise 

S 8h (0 =«%(+!) -<%(-l) 



1 , V k2 F+pi 

kd 



cosh — — 9(ef — eV)e 
kd 



-k F /k d 



kd 



kd 



V k d 



(74) 



Again we see the non-analyticity of this expression with respect to the driving field. In this uniformly-embedded wire 
model, expansion of the shot noise in powers of V is not valid at low fields. On the contrary, the shot noise becomes 
perturbative in 1/V at high fields, which could be termed the Schottky domain. The expression is perfectly calculable 
and there are two asymptotic cases of interest, 
(a) High fields, eV ^> £p: 



2el 1 



Pd 



-Pd/kd 



2el 1 



m*(l/ T ) 
eV~ 



(75) 



The high-field limit gives the full Schottky expression 2el with a correction, dependent on the wire parameters, which 
is asymptotically negligible. The same formal result also holds for any chosen value of / in the collisionless regime 
t — > oo, confirming that our model recovers the shot-noise behaviour of a monoenergetic flux. 

(b) Low fields, eV <C ep. The mean free path is A = rvp. Then (kp ±p 2 d )z — > kp[l± Ikd/Xkp — \(lkd/ Xkp) 2 ] and 



5 sh (0 



( I I 2 



(76) 



Our ID inelastic model again gives 2el at low fields, with no suppression in the ballistic limit A ^> I. There is, 
hovever, exponential decay of the shot noise as the wire length increases beyond the mean free path. The result is 
understandable as source and drain currents rapidly decorrelate with increasing i/A; if A oc r is made smaller while 
I is kept fixed, the exponential attenuation is broadly consistent with Monte Carlo results of Liu, Eastman, and 
Yamamoto.E3 Their more general simulation of ID low-field shot noise, which includes both elastic scattering and 
inelastic phonon emission, exhibits strong suppression in the inelastically dominated regime. 

In Fig. |^ we compare the ID shot noise normalised to 2el, for degenerate and classical conductors. Fig. ^(a) shows 
the results for a degenerate sample, as a function of current normalised to If — nevp. The plots are for a ran ge o f 
wire lengths from the ballistic limit I = 1CP 6 A, up to I — 10A. At low fields the intercepts at / = 0, given by Eq. (|7q), 
show their attenuation away from the ballistic limit. For higher currents, the shot noise quickly settles to the form 
given by Eq. (f75j). 

In Fig. ^(b) we plot, for comparison, the shot noise of a ID system whose carrier distribution is classical: A/^ q = 

oc exp(— ek/m*v 2 h ) where t>th = {ksT /m*)? . The current is normalised to / t h = ne^f^v^ and A = v^t.. Also 
plotted is the asymptotic form appearing in the first line of Eq. (|7^) . Again at low currents we see attenuation with 
increasing wire length, stronger than in Fig. |^(a). At higher currents there is the same rapid convergence to the 
asymptotic result (evident at surprisingly modest currents) as found in Fig. |^(a). 

From our comparison of Fermi-Dirac and Maxwell-Boltzmann versions of the model, we conclude that degeneracy 
contributes mainly at currents below Ip. For eV < ep the driving voltage cannot overcome the collective stability 
of the Fermi sea, and the zero-current correlations persist. For eV ^ £p appreciable redistribution of the particle 
occupancies suddenly becomes possible, with an initial dip in relative correlation strength. At higher fields most 
electrons move independently and ballistically, in the sense that r(u/)/n 3> I. The shot noise is then in the Schottky 
domain. 
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B. Two and Three Dimensions 



For higher dimensions we must include traces over the transverse degrees of freedom. Write A = 2R for v = 2 
and A = ttR 2 for v — 3 where R is the half- width (for a strip) or the radius (for a cylinder). As in the ID case, the 
condition / <C L x implies that the correlated BGF is well approximated by 

C k k'(q,q',w) = AS(q x - q' x )ip,y-i{q± - q' ± ; i?)C^,(q, w); 

refer to Appendix |c] for details. After integrating over the cross-sectional co-ordinates and applying Eq. ( |C9h| ) for 
C(°), Eq. © reads 



ne 2 A ( Ti\ 



2. kx 



k' x dk' x ^9(k x - k' x )e-^- k '^/ kd 
kd 



The shape factor F„(a), whose argument is a = hr\k x — k x \k±/(2m*kdR), has the form 
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which reduces to the expressions 



/d v q' 
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(77) 



(78a) 



(78b) 



F 3 (a) =6(1 -a) 



— arcsm a 
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(78c) 



In a wire of finite width, this function directly expresses the constraint on lateral motion of the carriers; it cross- 
couples, kincmatically, the transverse and longitudinal modes. Here there is none of the dynamical cro ss-co uplin g 
induced, for example, by elastic scattering. Kinematic suppression is inherent in the form of [Eqs. (C8) and (C9)], 
itself conditioned by the free-streaming operator in the Boltzmann equation. It is not surprising to find an entirely 
geometric source of shot-noise suppression in two and three dimensions (2D; 3D). This echoes, in part, Landauer's 
remarket on the need to sample more than just the longitudinal trajectories in any semiclassical c alcula tion. 

We process Eq. ( |77| ) along lines analogous to Eq. (|7^), going to the degenerate limit with Eq. (C14b) for Af and 
the fluctuation density An = m*kBTk F ~ 2 /fi 2 7r l/_1 . The first shot-noise component reduces to 
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(79a) 



wherep^ = (k 2 F — k 2 ^)? andpa, = (k 2 —p 2 d )^; the shape-factor arguments are a± = hT(k x ±p x )k±/(2m*kdR)- Similarly 
the second component is 
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(79b) 
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In Figs. H and [| we plot the shot noise in two and three dimensions for a range of wire geometries, as a function of 
current normalised to Ip — neAvp. The results in Fig. [5] are calculated for the thick-wire limit R 3> A; those in Fig. 
[| are for a thin wire, R — 0.3A. The same values of X/l are used as in Fig. Ei running down monotonically from the 
top. 

For both strips and cylinders in Fig. it is the behaviour of iS s h at higher values of / that comes to notice first: each 
curve merges with its asymptotic ID analogue illustrated in Fig. ||(b). As will shortly become clear, it is only in the 
thick-wire limit, and then only for high currents, that a ID treatment can in any sense mimic the exact calculations 
for higher dimensions [those conditions amount to setting F v sa 1 within Eq. (|79|)1. In fact, as one progresses from 
ID through 3D [Figs. ||(a), ||(a), and||(b)], the zero-current intercepts of the curves indexed by the same X/l undergo 
a marked and systematic increase in suppression. This shows (at least in the simple Drude model) that the ID 
calculation is a poor estimate of low-field noise in realistic geometries, even in the thick-wire limit. 

We come to the thin wires of Fig. 0. The uppermost curves, in the ballistic regime A > I, are unchanged from Figs. 
^| and ^. However, relative to Fig. g the remaining curves in Fig. ^ change dramatically as one moves further from 
the ballistic limit. There is now substantial shot-noise suppression throughout the whole range of /. For example, in 
Fig. ||(b) the longest 3D wire, I = 10A, has 6> s h/2e/ = 0.59 at the highest current / = lOIp, while its opposite number 
in Fig. 0(b) reaches the value 0.2; a threefold reduction. Calculations at much higher fields confirm the eventual 
recovery of full shot noise in keeping with Eq. (|75| ) . 

The outcome of spatially constrained carrier motion is the extensive suppression of shot noise over a wide range 
of the current. We stress that the effect is implicit in the generic structure of the Boltzmann equation, and its 
propagators, for 2D and 3D [see Eq. (C8)]; it is simply absent in ID. Therefore, this mode of suppression cannot be 
simulated by any one-dimensional scheme. 

In a 3D wire the Fermi wavelength is 0.05 nm at metallic electron densities. If A — 50 nm, typical for strong inelastic 
scattering, a wire of width ~ 30 nm would exhibit kinematic shot-noise suppression at large currents, providing that it 
was not masked by collisional Coulomb effects. In future we plan to assess the latter quantitatively; the comparative 
action of nonequilibcium screening, in 2D versus 3D, should itself be an interesting window on how dimensionality 
affects fluctuations rl 



V. SUMMARY 



We have described and implemented a nonperturbative microscopic formalism for current fluctuations in metallic 
systems, down to the mesoscopic scale, within the ambit of semiclassical theory. Our strategy for incorporating 
fermion correlations into the Boltzmann picture safeguards the conservation laws at both the single-particle level 
and at the level of dynamic two-particle processes, the key to nonequilibrium current noise. In particular we have 
derived the nonlinear analogue of the fluctuation-dissipation theorem. It should also be straightforward to include 
semiclassical electron-electron scattering in our description of transport and noise. 

Our formalism's calculability stems from the flexible structure of the Boltzmann-Green functions. These serve as 
semiclassical propagators of the electronic Fermi-liquid correlations, mapping them uniquely to the correlations of the 
nonequilibrium system. We have demonstrated their usefulness in shedding light on the physics of high-current shot 
noise, and on the importance of treating dimensionality correctly in constricted mesoscopic samples. 

The present account of nonequilibrium fluctuations raises a variety of interesting questions, practical and abstract. 
By far the most salient is the relation between thermal noise and shot noise; our claim that the two are thermodynam- 
ically incommensurate should be easy to test. Thermal-noise measurements on a gated two-dimensional mesoscopic 
wire, defined on a III-V heterojunction, ought to give a strong gate-voltage-dependent signature of Thomas-Fermi 
suppression from self-confinement of the carriers in their quantum well. Measurements of the shot noise in the same 
structure should give no such signature. 

We end with just two out of many theoretical issues. The first is the role of non-analyticity of the BGF solutions 
at low f)*jdS|, and the implications for semiclassical linear response. The phenomenon is well known in uniform 
systems J29'c3 where the free-streaming part of the Boltzmann operator is manifestly anomalous in its vanishing with 
the applied field. While we have no firm information on whether the BTE for nonuniform systems shares this 
behaviour, we point out that our shot-noise results, obtained in a spatially inhomogeneous model (albeit weakly so), 
are certainly non-analytic in the applied voltage. The low-field asymptotics of the general Boltzmann equation have 
an obvious bearing on how semiclassical noise is to be calculated, and their clarification would be a significant advance. 

Second, there is the status of the_BGF approach within quantum kinetics. In terms of, say, the Kadanoff-Baym 
analysis of the Boltzmann equationjfj our semiclassical equation of motion for the fluctuations should emerge from 
the quantum evolution of the particle-hole amplitudes in the long- wavelength limit ,Cj much as the ordinary BTE is 
distilled from the long-wavelength dynamics of the density matrix for a single particle. 
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APPENDIX A: TYPES OF SCREENING IN NONEQUILIBRIUM SYSTEMS 



We begin this Appendix with some basic properties of collision-mediated Coulomb screening, which influences both 
thermal and shot noise. We end with Thomas-Fermi screening, which influences thermal noise alone. Keeping in 
mind the high-field application, we base our calculations on the collision time model (see the following Appendices). 
Our results are specific to a bulk system with inelastic scattering; we stress both their illustrative intent and their 
inapplicability to the elastic diffusive regime, where the Boltzmann propagators and the collisional screening effects 
are very different. 

We take the Fourier transform of the collision-mediated Coulomb operator [see Eq. (f43|)] using the form of the 
correlated BGF in the bulk limit, Eq. ( p!2| ). The (spin-independent) function T is 



Tkk'(q; w) = fW k k' 



2e 



a/k» 

dk" 



•^c(q)5Z rk '" k '( c i!' J )' 



(Al) 



where the electron field is eScil) — ~*gVc(<?) and Vc(q) is the Coulomb potential transform. In three dimensions 
there is a complication owing to the long-range Coulomb tail; by assumption, all fields are shorted out beyond the 
system boundaries. We model this constraint by introducing a cutoff, so that Vc(<?) = ^e 2 /e(q 2 + K 2 ) where k~ 1 > f2s 
represents the characteristic length scale beyond which the fields are zero. 

The trace j co \\(q,uj) = (r(q, w;k')) is independe nt o f wave vector k'. This follows from the decoupling of the 
internal summations over kinematic variables in Eq. ( |Al[ ) , since the Coulomb field depends only on q. Summing both 
sides of Eq. (jAl]) we get 

2e x-^ v-v , x dfk 



7cou(q,o;) = 1- jjj2^^C^(q,w)-^-£ c (q)7coii(q,w) 
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dk' 
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dfv 
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(A2) 



There is now a closed form for the Coulomb operator: 

.87re 2 7 co n(q,w) 



Lkk'(q, w) = fWkk' + i- 



E c S-(q,-)q 



dfk» 
dk" ' 



(A3) 



suggesting a possible approximation for the convolution of T with a typical distribution F, namely 

^^r kk /(q,ixi)Fk' ~ 7coii(q,^)i ; k . 



(A4) 



This is exact in the q — > limit and also reproduces the exact relation {(T(q,ui)F'))' — 7coll(q> u))(F). The Ansatz, 

which amounts to the decoupling C k fc ^//(C^^(k')) ~ Fk/(F), tends to wash out the sharp features of the integrand in 
Eq. ©. 

We evaluate Eq. (A2) for the Drude model in the zero- field limit. Using Eq. ( |C13| ), the trace of the correlated 
BGF over its leading wave vector is 
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There is no contribution to Eq. (|A2|) from the second right-hand term of Eq. (A5), and so 
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7coii(g,0) 



1 + gTF ( x _ arctan(Ag) \ 
q 2 + k 2 \ Xq J 



(A6) 



where the Thomas-Fermi wave vector is defined by q^ F — kite 2 An/ tk B T = 4e 2 m*A;F /"nth 2 , and A = rvp is the mean 
free path. In the small-g limit the suppression factor becomes 



7coii(g,0) 



k 2 + (1 + AnaT/e)q 2 



(A7) 



(er is the conductivity), showing that there is no long-range suppression from collision- mediated screening if the 
asymptotic state of the leads pins the electric fields to zero there. Put differently, microscopic scattering preserves 
global neutrality; Thomas-Fermi screening may not, since it redistributes charge. This calls for the inclusion of buffer 
zones at the system boundaries (in practice, several units of <?tf) t° ensure that the fields beyond the system remain 
evanescent. 

Mesoscopically, collisional screening should be significant. In a metal ~/ co \\(q,0) can be very small; for example, in 
silver at 77 K its minimum is roughly 10~ 7 and its value does not rise to one half until q^ 1 = q^ F = 0.06 nm, that 
is, far below any mean free path and out of the semiclassical domain. In heavily doped GaAs with carrier density 
n = 10 18 cm -3 , comparable figures are 7 co n = 0.01 for the minimum and 7 co ii = 0.5 at q^ 1 — 10 nm. While these are 
guideline figures for a simple model in its approach to equilibrium, they hint at a strong role for collisional screening 
suppression in high-field shot noise. 

We return to Thomas-Fermi screening, associated with the contact potentials. This is a primary source of thermal- 
noise suppression, a thermodynamic effect free of any collision processes. We outline its behaviour in a bulk jellium 
conductor contacted by leads made of different jellium. The combined system is treated as an electron gas closed with 
respect to carrier exchange, and satisfying periodic boundary conditions. □ We must also take explicit account of the 
reservoir's electrostatic potential u r . 

If the subsystems are macroscopic, a term such as the second one in the right-hand-side sum of Eq. (|4l]) goes to 
its asymptotic mean u c ; the third term is negligible to the same order. Including the explicit offset from it r , Eq. ( |39| ) 
generalises to 
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(A8) 



where all primed quantities refer to the reservoir. A complementary relation holds for 7^(5). Consider the conductor. 
In the long-wavelength limit we have 

Vc(9)^ = %■ 
k B T q l 



Putting this into Eq. (A8) together with its reservoir counterpart, and taking q <C <Ztf,9tf: we obtain the coupled 
equations 



u r An 



9tf 



k B Tn 

giving the limiting solutions 



70(g) - Utf + |^« 2 Vo(g) = ? 2 ; 7c(g) + i c [q) 



(A9) 



7C = 1 _ itf q tf and ic = 1 , 5tf g TF 



gxF + q'rF 



9tf + 9'tf 



(A10) 



In this system, a conductor that is more metallic than the reservoir has ^tf > 9tfi leading to 7c < 1 < 7^7 < 2. Its 
thermal noise (were it accessible from outside) would thus undergo suppression. Conversely, a relatively less metallic 
sample would display enhanced thermal noise. If (?tf "C q' TF , as in a lightly doped bulk semiconductor in contact with 
a metal, then 7c approaches its maximum of two. 
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A closed model provides an unrealistic picture of actual transport; it may be indicative, but by no means definitive. 
In particular, for truly open reservoirs, the potential u r cannot be accessed in isolation from the chemical potential 
as a whole. This means that the reciprocity between jc and 7^ , characteristic of the closed model, does not hold. 
Instead, the identity j' c = 1 is a necessary boundary constraint. 

For an open mesoscopic conductor, details of the contact-potential effects on thermal noise will be sensitive both to 
the topology of its boundary conditions, and to its internal electronic structure. The action of suppression (or indeed 
of enhancement) is a problem for further study. 



APPENDIX B: UNIFORM DRUDE MODEL 



We derive the dynamical fluctuation structure for a single parabolic conduction band with uniform electron density 
n and constant mobility fi e = er/m*, where r is the spin- independent collision time and m* the effective mass. The 
system is driven by a uniform field E = E = — Eit acting in the negative (drain to source) direction. We take variations 
which are homogeneous over the sample region, so that the fluctuations of interest have no spatial dependence. 

The Boltzmann equation in the model is 



d_ 

dt 



eE d 
h dk x 



hit) 
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eq 
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(BI) 



Since the Boltzmann operator is linear, the fluctuation structure is qualitatively similar to that for elastic scattering 
[differences arise from the inhomogeneous term in / eq , notably in the behaviours of R(t) and A/(t)]. We solve Eq. 
( [BID by Fourier transforms in reciprocal space, so that the transform F p = f2 _1 ^ k /k exp(ik-p) of the steady-state 
distribution takes the form 



F 



F — 
p ~ F n eq 1 



P 



ikdPx 



(B2) 



where kd = eEr/fi and Fq = \ (f) per spin state. Note that, while a formal distinction is made between Fq and F^, 
the physical normalisation is always Fq = F, 



eq 



The transform of the dynamic BGF, TZ pp >(u>) = ft 2 J2u k'R-kk'(^) exp[i(k-p — k'-p')], has the equation 



-iuT - ik d p x + 1] Hp p i (uj) =t8(p- p') + 
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For p = this leads to 



Ha P '{u) = - 



S(P') 



i(uj + irj) 



(B3) 



(B4) 



On the other hand, the uncorrelated component of 7Z PP > scales with the steady-state solution F p [in a collision-time 
model the asymptotic form F p /^n replaces AF p /^An]. Denoting the correlated part by C pp i and recalling that the 
uncorrelated part exhausts the normalisation of 1Zq p > , we obtain 



Ti ppl {u)) = C pp ,{uj) 



S(P') F p 

i(ld + IT]) Fq 



(B5) 



When the above is put together with Eqs. (B2)-(B4) we arrive, after some algebra, at the explicit formula for the 
correlated propagator: 



6(p-p')-^8(p>) 

f - ikdPx - iut 



(B6) 



We can use Eq. (B6) directly to evaluate both dissipative and non-dissipative contributions to the noise. Using the 
reciprocal-space representation v «-* —i(H/m*)d/dp, the power density P of Eq. (|60| ) is 
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aE 2 



(B7) 



The Drude conductivity a — nep e appears when we apply the relation {— d 2 AF p iq /dp 2 .} Q = (k 2 Af eq ) 
m*kBTn/2h 2 to the middle line of the equation. A contributio n co ntaining (v x Af cq ) = vanishes trivially. 



The hot-electron spectral density S g in the static limit [Eq. ( J64a| )] is calculated similarly: 

feE-x) 3 



p^O 



3 E 3 r 2 h \ 


' d 


1 


( d -ipx 


m* 2 ) 


dp x 


1 - ikdPx \dp x 1 - i 


' d F p /F 




' d -ip' x AFp 


dp x 1 - i 


kdPx_ 


p^O 


dp' x 1 - ik d p' x 



p^O 



(B8) 



We evaluate this with the help of the relations AFp q = |An and {dF p /dp x }p->o — ikdFo, the latter following from 
Eq. @. The result is 



,-- * 2 77*4 ' 

b g = am [i e hj 
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APPENDIX C: WEAKLY NONUNIFORM DRUDE MODEL 



We derive the spatio-temporal correlations within the Drude model of the preceding Appendix. The problem is to 
calculate the propagation of a single electron added to N uniformly distributed electrons at a specific point in the 
sample at t — 0. This constitutes a weak inhomogeneity. 

Since the scattering is spin-independent, we consider a zero-spin model with effective density ^n. The equation of 
motion for the dynamical propagator in the frequency domain is 



d , d 1 . 
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or ak T 
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Now define the Fourier transform of the propagator, 

R kk /(q,q» = /Vr /"dV R QCe 'H exp[-i(q-r - q'-r')]. 
Jo Jn 



Eq. (CI) becomes 
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where 



(C4) 



For a three-dimensional system with cylindrical symmetry about the x-axis we write q = (q x ,q±) and f2 = irR 2 L x 
where L x is the sample length and R its radius. The function ip^ can be written as 



¥>3(q) = ¥>ite; |ix)¥'2(q_L;-R) = 



\L x q x hRqx 
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with Ji(u) the first-order Bessel function. For a two-dimensional strip of width 2R, the role of ips is assumed by the 
product (pi{q x ;^L x )ipi(q±;R). 

The correlated dynamical propagator C associated with R has the form 



C kk ' (q, q» = R kk ' (q, q', u) - (R(0, q', w; k'))w(q) { 



fa 



= R k k'(q,q',w) 
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wherein-^R) is evaluated by summing over k on both sides of Eq. ( |C3| ) in the limit g — ► 0. The equation of motion 
for C isN 
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Terms ~ singular in the static limit, cancel identically by the fact that / is the Bqllpxiann solution for the 

uniform steady state. Equation ( |C7j ) is solved at arbitrary fields with the integrating factorEj 
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in the low-field limit, the non- analytic character of X is clear from the occurrence of 1/kd oc 1/E in its exponent. 
Using X we first generate the auxiliary propagator C^ ^ satisfying 
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The expression for C^ ^ is 
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which furnishes the complete solution to Eq. (C7) as 
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For our study of shot noise it is convenient to consider the propagator integrated over q', namely 

C&(q,") = / |^;C kk ,(q,q>), 



(Cll) 



equivalent to the bulk solution in an infinitely wide conductor, R — > oo. The form of C^ follows from Eq. (CIO). It 
is 
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In the zero-field limit this goes to 
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req 
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Finally, for reference, we also record the expression for the uniform distribution A/k, needed in the shot-noise 
application: 
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In the degenerate limit this becomes 
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where pj_ = — £;5_) 2 . In one dimension, set k± = 0. 



2 



5 



t Electronic address: Fred.Green@tip.csiro.au 
' Electronic address: mukunda.das@anu.edu.au 
1 We refer to calculations based on the microscopies of degenerate systems rather than predominantly numerical methods such 
as Monte Carlo; see for example P. Tadyszak, F. Danneville, A. Cappy, L. Reggiani, L. Varani, and L. Rota, Appl. Phys. 
Lett. 69, 1450 (1996). 

D. Pines and P. Nozieres, The Theory of Quantum Liquids (Benjamin, New York, 1966). 

3 N. W. Ashcroft and N. D. Mermin, Solid State Physics (Holt, Rinehart, and Winston, New York, 1976). 

4 M. Reznikov, M. Heiblum, H. Shtrikman, and D. Mahalu, Phys. Rev. Lett. 75, 3340 (1995). 
A. Kumar, L. Saminadayar, D. C. Glattli, Y. Jin, and B. Etienne, Phys. Rev. Lett. 76, 2778 (1996). 

6 A. H. Steinbach, J. M. Martinis, and M. H. Devoret, Phys. Rev. Lett. 76, 3806 (1996). 

7 R. J. Schoelkopf, P. J. Burke, A. A. Kozhevnikov, D. E. Prober, and M. J. Rooks, Phys. Rev. Lett. 78, 3370 (1997). 

8 H. Pothier, S. Gueron, N. O. Birge, D. Esteve, and M. H. Devoret, Phys. Rev. Lett. 79, 3490 (1997). 

9 R. de Picciotto, M. Reznikov, M. Heiblum, V. Umansky, G. Bunin, and D. Mahalu, Nature 389, 162 (1997). 

10 L. Saminadayar, D. C. Glattli, Y. Jin, and B. Etienne, Phys. Rev. Lett. 79, 2526 (1997). 

11 For a comprehensive overview, see M. J. M. de Jong and C. W. J. Beenakker in Mesoscopic Electron Transport, edited by 
L. P. Kouwenhoven, G. Schon, and L. L. Sohn, NATO ASI Series E (Kluwer Academic, Dordrecht, 1997). 

12 G. B. Lesovik, JETP Lett. 66, 1243 (1989). 

13 M. Biittiker, Phys. Rev. Lett. 65, 2901 (1990); Phys. Rev. B 46, 12485 (1992). 

14 C. W. J. Beenakker and M. Biittiker, Phys. Rev. B 46, 189 (1992). 

15 Th. Martin and R. Landauer, Phys. Rev. B 45, 1742 (1992). 

16 K. E. Nagaev, Phys. Lett. A 169, 103 (1992); Phys. Rev. B 52, 4740 (1994). 

17 M. J. M. de Jong and C. W. J. Beenakker, Phys. Rev. B 51 16867 (1995). 

18 R. Landauer, IBM J. Res. Dev. 1, 233 (1957). 

19 Sh. M. Kogan and A. Ya. Shul'man, Zh. Eksp. Teor. Fiz. 56, 862 (1969) [Sov. Phys. JETP 29, 467 (1969)]. 

20 The transition from coherent to incoherent prescriptions for noise has been studied by A. Shimizu and M. Ueda, Phys. Rev. 
Lett. 69, 1403 (1992). 

21 R. Landauer, Phys. Rev. B 47, 16427 (1993). 

22 R. C. Liu, P. Eastman, and Y. Yamamoto, Solid State Commun. 102, 785 (1997). 

23 R. Landauer, Physica (Amsterdam) 227B, 156 (1996). 

24 D. K. Ferry and S. M. Goodnick, Transport in Nanostructures (Cambridge University Press, Cambridge, U. K., 1997). 

25 C. J. Stanton and J. W. Wilkins, Phys. Rev. B 35, 9722 (1987); Phys. Rev. B 36, 1686 (1987). 

26 B. L. Altshuler, L. S. Levitov, and A. Yu. Yakovets, Pis'ma Zh. Eksp. Teor. Fiz. 59, 821 (1994) [JETP Lett. 59, 857 (1994)]. 

27 T. Gonzalez, C. Gonzalez, J. Mateos, D. Pardo, L. Reggiani, O. M. Bulashenko, and J. M. Rubf, Phys. Rev. Lett. 80, 2901 
(1998). 

28 G. Rickayzen, Green's Functions and Condensed Matter (Academic Press, London, 1980). 

29 S. V. Gantsevich, V. L. Gurevich, and R. Katilius, Riv. Nuovo Cimento 2, 1 (1979). 

30 C. J. Stanton, Ph.D. thesis, 1986 (Cornell University, unpublished). 



23 



31 F. Green and M. J. Chivers, Phys. Rev. B 54, 5791 (1996). 

32 Y. Naveh, D. V. Averin, and K. K. Likharev, Phys. Rev. Lett. 79, 3482 (1997). 

33 T. Gonzalez, O. M. Bulashenko, J. Mateos, D. Pardo, and L. Reggiani, Phys. Rev. B 56, 6424 (1997); O. M. Bulashenko, J. 
Mateos, D. Pardo, T. Gonzalez, L. Reggiani, and J. M. Rubf, Phys. Rev. B 57, 1366 (1998). 

34 K. E. Nagaev, Phys. Rev. B 57, 4628 (1998). 

35 C. M. Van Vliet, IEEE Trans. Electron Devices 41, 1902 (1994). 

36 J. P. Nougier, in Physics of Nonlinear Transport in Semiconductors, edited by D. K. Ferry, J. R. Barker, and C. Jacoboni 
(Plenum, New York, 1980), p. 415 ff.; see also J. P. Nougier, IEEE Trans. Electron Devices 41, 2034 (1994). 

37 L. Reggiani, P. Lugli, and V. Mitin, Phys. Rev. Lett. 60, 736 (1988). 

38 The Poisson equation is always defined in three dimensions. To interpret Eq. (^) appropriately when v < 3, the electron 
density (/) must be understood to carry a (separable) factor in the 3 — v transverse space co-ordinates. Thus, for a transport 
problem confined strictly to two dimensions, the Poisson source term contains (/(r)) = 8(z){f(r±)), where z is orthogonal 
to the plane (r±;z — 0). On the other hapd, the stabilising background distribution n + (r) can be fully three-dimensional, 
as in a modulation-doped heterostructureO 

39 F. Green, Phys. Rev. B 54, 4394 (1996). 

40 D. K. C. MacDonald, No ise and Fluctuation s: an Introduction (Wiley, New York, 1962). 

41 



F. Green and M. P. Das, cond-mat/9709142 



' P. M. Bakshi and E. P. Gross, Ann. Phys. 24, 419 (1968). 

' R. Kubo, M. Toda, and N. Hashitsume, Statistical Physics II: Nonequilibriurn Statistical Mechanics (2nd ed., Springer, 
Berlin, 1991), pp 199 and 200. 

The BGF description does produce a smooth crossover from thermal to shot-noise- like behaviour in some regimes, but always 
scaled by the degeneracy factor An/n. See our Ref. £y|; for the classical limit see Stanton, Ref. 

' L. P. Kadanoff and G. Baym, Quantum Statistical Mechanics (Benjamin, Reading, MA, 1962). 

' Sh. M. Kogan, Phys. Rev. A 44, 8072 (1991). 



24 



FIG. 1. Zero-frequency spectral density of nonequilibrium thermal noise in a uniform, two-dimensional electron gas in 
GaAs, plotted as a function of the external field from T — to 900 K in 150 K steps. Normalisation is to the equilibrium 
Johnson- Nyquist value. At low temperature, degeneracy sets the scale of the contribution from nonequilibrium electron heating. 
At high temperature, the hot-electron component shifts up towards higher fields as the equilibrium component gains dominance. 
Dot-dashed line: thermal noise at 300 K. 

FIG. 2. (a) Shot noise in the Drude model of a degenerate one-dimensional wire, as a function of current in units of the 
Fermi current and normalised to full shot noise. Each curve is for a fixed ratio X/l of mean free path to length of wire; the 
same set of ratios is used for all subsequent figures. At high currents and in the ballistic limit (upper curves), the shot noise 
tends to its full value. At low currents and away from the ballistic limit, degeneracy inhibits the natural tendency of the shot 
noise to exponential suppression with increased sample length, (b) Full line: shot noise of classical carriers, as a function of 
current in units of the thermal current. The attenuation at low currents is much more pronounced than in (a). Dots: high-field 
asymptote defined by Eq. ((r5]). Note how both classical and fermionic results rapidly assume this form at higher currents. 

FIG. 3. Shot noise in the degenerate Drude model for (a) very wide two-dimensional strips, and (b) three-dimensional 
cylinders of very large radius. The behaviour at low currents differs significantly from Fig. ^|(a), with progressively greater 
suppression at longer wire lengths and higher dimensionality. At high currents the asymptotic behaviour is identical with that 
in one dimension. 

FIG. 4. Shot noise in the degenerate Drude model for (a) narrow strips, and (b) thin cylinders. Away from the ballistic 
limit (topmost curves) there is remarkable shot-noise suppression over the entire range of the current. This effect is inherent in 
the kinematic term of the higher-dimensional Boltzmann equation; it cannot be simulated by a one-dimensional approximation. 
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